A fast and efficient algorithm for Slater determinant updates in quantum Monte Carlo simulations.
نویسندگان
چکیده
We present an efficient low-rank updating algorithm for updating the trial wave functions used in quantum Monte Carlo (QMC) simulations. The algorithm is based on low-rank updating of the Slater determinants. In particular, the computational complexity of the algorithm is O(kN) during the kth step compared to traditional algorithms that require O(N(2)) computations, where N is the system size. For single determinant trial wave functions the new algorithm is faster than the traditional O(N(2)) Sherman-Morrison algorithm for up to O(N) updates. For multideterminant configuration-interaction-type trial wave functions of M+1 determinants, the new algorithm is significantly more efficient, saving both O(MN(2)) work and O(MN(2)) storage. The algorithm enables more accurate and significantly more efficient QMC calculations using configuration-interaction-type wave functions.
منابع مشابه
A sparse algorithm for the evaluation of the local energy in quantum Monte Carlo
A new algorithm is presented for the sparse representation and evaluation of Slater determinants in the quantum Monte Carlo (QMC) method. The approach, combined with the use of localized orbitals in a Slater-type orbital basis set, significantly extends the size molecule that can be treated with the QMC method. Application of the algorithm to systems containing up to 390 electrons confirms that...
متن کاملSufficient condition for absence of the sign problem in the fermionic quantum Monte Carlo algorithm
Quantum Monte Carlo sQMCd simulations involving fermions have a notorious sign problem. Some wellknown exceptions to the auxiliary field QMC algorithm rely on the factorizibility of the fermion determinant. Recently, a fermionic QMC algorithm fC. Wu, J. Hu, and S. Zhang, Phys. Rev. Lett. 91, 186402 s2003dg has been found in which the fermion determinant may not necessarily be factorizable, but ...
متن کاملKinetic Monte Carlo Simulation of Strained Heteroepitaxial Growth with Intermixing
An efficient method for the simulation of strained heteroepitaxial growth with intermixing using kinetic Monte Carlo is presented. The model used is based on a solid-on-solid bond counting formulation in which elastic effects are incorporated using a ball and spring model. While idealized, this model nevertheless captures many aspects of heteroepitaxial growth, including nucleation, surface dif...
متن کاملAccuracy of electronic wave functions in quantum Monte Carlo: The effect of high-order correlations
Compact and accurate wave functions can be constructed by quantum Monte Carlo methods. Typically, these wave functions consist of a sum of a small number of Slater determinants multiplied by a Jastrow factor. In this paper we study the importance of including high-order, nucleus-three-electron correlations in the Jastrow factor. An efficient algorithm based on the theory of invariants is used t...
متن کاملMonte Carlo Simulations for Quantum Field Theories InvoMng Fermions
We present a new variant of a Monte Carlo procedure for euclidean quantum field theories with fermions. On a lattice every term contributing to the expansion of the fermion determinant is interpreted as a configuration of self-avoiding oriented closed loops which represent the fermionic vacuum fluctuations. These loops are related to Symanzik's polymer description of euclidean quantum field the...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- The Journal of chemical physics
دوره 130 20 شماره
صفحات -
تاریخ انتشار 2009